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Abstract 

Background: Breast cancer survivors, particularly those treated with chemotherapy, are at significantly increased risk 
for long-term cognitive and neurobiologic impairments. These deficits tend to involve skills that are subserved by 
distributed brain networks. Additionally, neuroimaging studies have shown a diffuse pattern of brain structure 
changes in chemotherapy-treated breast cancer survivors that might impact large-scale brain networks. 

Methods: We therefore applied graph theoretical analysis to compare the gray matter structural networks of 
female breast cancer survivors with a history of chemotherapy treatment and healthy age and education matched 
female controls. 

Results: Results revealed reduced clustering coefficient and small-world index in the brain network of the breast 
cancer patients across a range of network densities. In addition, the network of the breast cancer group had less 
highly interactive nodes and reduced degree/centra I ity in the frontotemporal regions compared to controls, which 
may help explain the common impairments of memory and executive functioning among these patients. 

Conclusions: These results suggest that breast cancer and chemotherapy may decrease regional connectivity as 
well as global network organization and integration, reducing efficiency of the network. To our knowledge, this is 
the first report of altered large-scale brain networks associated with breast cancer and chemotherapy. 



Background 

Breast cancer is one of the most frequently diagnosed 
cancers and the leading cause of cancer death among 
females, accounting for 23% of the total cancer cases [1]. 
In the last decade, several neuropsychological studies 
have shown the negative influence of breast cancer and 
chemotherapy on various cognitive skills with executive 
function and memory impairments being the most com- 
mon [2,3]. These deficits have been reported both prior 
and following chemotherapy with evidence showing 
increased and/or more severe cognitive changes in 
breast cancer patients treated with chemotherapy [4-9]. 
Neuroimaging studies corroborate these findings by 
showing changes in both brain structure and function 
associated with breast cancer and chemotherapy [10]. 
However, it is currently unknown whether breast cancer 
and chemotherapy affect large-scale brain networks. 
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There are several lines of evidence suggesting that 
breast cancer may negatively impact whole brain net- 
works. First, while the mechanisms by which breast can- 
cer and its treatments affect cognitive function are 
largely unknown, possible candidates include neurotoxic 
effects of chemotherapy, oxidative damage, cytokine dys- 
regulation and individual variation in genes related to 
neural repair and/or plasticity [2,3,11,12]. These candi- 
date mechanisms are likely to have diffuse effects on 
brain structure. Second, neuroimaging studies indicate 
that breast cancer survivors show altered brain structure, 
which would disrupt large-scale networks [13]. Specific- 
ally, these patients demonstrate reduced gray matter in 
bilateral frontal, temporal, cerebellar, thalamic, and cin- 
gulate regions as well as decreased white matter integrity 
in corpus callosum, frontal, and temporal white matter 
tracts [14-17]. Third, the specific cognitive domains that 
tend to be most commonly impaired in breast cancer in- 
volve executive functions and memory, as noted above. 
These skills are known to be subserved by distributed, 
integrated neural networks [18]. Finally, cognitive im- 
pairment following breast cancer often tends to be quite 
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subtle [2], potentially suggesting a more diffuse brain in- 
jury. In the current study, we investigated whether 
breast cancer and chemotherapy are associated with 
alterations in large-scale structural brain networks. 

Recent graph-theoretical analyses have consistently 
shown that brain structural networks in healthy indivi- 
duals have small- world characteristics [19]; an architec- 
ture that has dense local clustering of connections 
between neighboring nodes with short path length be- 
tween any pair of nodes due to the existence of relatively 
few long-range connections [20]. These characteristics, 
shared by various biological systems, reflect a network 
that is simultaneously highly segregated and integrated 
and allows for higher, more efficient rates of information 
processing and learning than random networks [21]. 

Since small-world characteristics were described quan- 
titatively for brain structural networks, there have been 
multiple graph-theoretical studies seeking to assess the 
structural correlation networks constructed from regional 
gray matter volume, cortical thickness and surface area 
[12,20,22,23]. The unique feature of graph-theoretical 
analysis, compared with the more traditional univariate 
neuroimaging approaches, is that it can directly test the 
differences in topological parameters of the brain net- 
work such as small-worldness, highly connected hubs 
and regional network parameters. Whereas univariate 
neuroimaging approaches have typically shown limited 
correlations with cognitive function and dysfunction, net- 
work parameters may provide a more robust model of 
cognitive status [24,25]. Recent graph- theoretical studies 
have illustrated an alteration of arrangements in struc- 
tural correlation networks associated with normal aging, 
multiple sclerosis, Alzheimer's disease, schizophrenia and 
epilepsy [26-30]. 

In the present study, we applied graph theoretical ana- 
lyses to compare magnetic resonance imaging (MRI)- 
based gray matter correlation networks of female breast 
cancer patients treated with chemotherapy and female 
healthy controls. Considering the lines of evidence 
regarding a diffuse pattern of gray matter atrophy in 
breast cancer patients, we hypothesized that such altera- 
tions should be reflected in small-world characteristics 
of the brain structural correlation network. We also 



examined the between group differences in highly con- 
nected hubs as well as in regional network measures 
such as node betweenness and degree. 

Methods 

Participants 

Thirty-seven female primary breast cancer survivors 
(BC) age 43-67 years who were treated with systemic 
chemotherapy (mean time since treatment = 4.5 ± 3.4 
years) and 38 age-matched healthy controls (CON) par- 
ticipated in the study. BC survivors were recruited via 
the Army of Women (http://www.armyofwomen.org/), 
community-based BC support groups and local media 
advertisements. Healthy controls were recruited via the 
Army of Women and local media advertisements. 
Groups were matched for age, education, premorbid 
cognitive functioning and minority status although there 
were significantly more post-menopausal women in the 
BC group (Table 1). Premorbid cognitive functioning 
was measured using the Information subtest of the 
Wechsler Adult Intelligence Scale 4th Edition. This subt- 
est measures fund of general knowledge and is believed 
to be influenced by previous experience and quality of 
education, making it a good estimate of premorbid abil- 
ity [31]. Participants in the BC group were free from dis- 
ease and had no history of relapse or recurrence at the 
time of evaluation. Chemotherapy regimens included 
AC (Doxorubicin, Cyclophosphamide) (N = 9), ACT 
(Doxorubicin, Cyclophosphamide, Paclitaxel) (N = 16), 
ACF (Doxorubicin, Cyclophosphamide, 5-Fluorouracil) 
(N = 2), CT (Cyclophosphamide, Paclitaxel) (N = 6), CMF 
(Cyclophosphamide, Methotrexate, 5-Fluorouracil) 
(N = 2), CTF (Cyclophosphamide, Paclitaxel, 5-Fluorour- 
acil) (N = l) and AC + CMF (N=l). Additionally, 15 
women were treated with tamoxifen and 25 received ra- 
diation therapy. None of the BC patients had cerebral 
metastasis or had been treated for that. Participants were 
excluded for neurologic, psychiatric or medical condi- 
tions known to affect cognitive function (e.g. learning 
disability, traumatic brain injury, chronic depression) as 
well as any MRI contraindications (e.g. non-MRI safe 
implants). This study was approved by the Stanford 



Table 1 Demographic data for the breast cancer and healthy control groups 





BC (N = 37) 


CON (N = 38) 


t or Chi Sq. 


p-value 


age 


54.2 (6.1) 


55.5 (9.0) 


.737 


.46 


education 


16 (2.8) 


17(2.6) 


1.61 


.11 


minority status 


8.6% 


10.5% 


2.27 


.69 


post menopause 


86.7% 


55.9% 


7.24 


.007 


premorbid cognitive status* 


13.0 (6.0) 


13.8 (2.7) 


.774 


.44 



^measured using the Information subtest of the Wechsler Adult Intelligence Scale 4th Edition. 
Data are shown as mean (standard deviation) except where noted. 
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University Institutional Review Board and all partici- 
pants provided informed consent. 

MRI data acquisition and preprocessing 

MRI scanning was performed on a GE Discovery MR750 
3.0 Tesla whole body scanner (GE Medical Systems, Mil- 
waukee, WI). High- resolution Tl -weighted images were 
acquired with 3D spoiled gradient echo pulse sequence 
using the following parameters: TR = 8.5 ms, TE = 3.396, 
TI = 400 ms, flip angle = 15, FOV = 220 mm, number of 
excitation =1, acquisition matrix = 256 x 192. Totally, 
124 contiguous coronal slices were with in-plane reso- 
lution of 0.859 mm x 0.859 mm. 

Image preprocessing was performed using Statistical 
Parametric Mapping 8 (SPM8; Wellcome Department of 
Cognitive Neurology, London, UK; http://www.fil.ion. 
ucl.ac.uk/spm/). The images were initially segmented 
into gray matter (GM), white matter, and cerebrospinal 
fluid images based on the ICBM Tissue Probabilistic 
Maps (http://www.loni.ucla.edu/ICBM/ICBM Tissue- 
Prob.html). A study-specific a priori probability map of 
GM was created from the modulated spatially normal- 
ized segmented GM images using the Template-O-Matic 
(TOM8) toolbox [32]. Next, the custom priors were af- 
fine-registered to the standard Montreal Neurological 
Institute (MNI) space. 

Voxel-based morphometry (VBM) preprocessing steps 
were performed using the VBM8 toolbox for SPM8 
(http://dbm.neuro.uni-jena.de/vbm/) and this involved 
segmentation of MR images into GM tissue segments 
using the standard unified segmented and non-linearly 
warping of the tissue segments to the GM study- specific 
customized template model [33]. Images were then 
modulated to ensure that relative volumes of GM were 
preserved following the spatial normalization procedure. 
Sample homogeneity was checked to identify any out- 
liers in the study population. Data for one participant in 
the BC group were excluded from the analysis because 
of a covariance below 2 SD and confirmation of visual 
motion artifact. 

ROI extraction 

There are different nodal definition methods in brain 
network analysis. While the results might be affected by 
the choice of panellation scheme, recent evidence 
showed that the results of between-group comparison 
remain intact regardless of the applied panellation 
scheme [34]. We generated 90 cortical and subcortical 
regions of interest (ROIs), excluding the cerebellum, 
from the Automated Anatomical Labeling (AAL) atlas 
using the WFU PickAtlas Toolbox [35]. The ROIs were 
identical to those used in a previous graph analysis study 
of typical brain development by Fan and colleagues [22]. 
These AAL ROIs were resliced to the same dimension 



as that of tissue segmented images obtained from the 
VBM preprocessing step. The ROIs were subsequently 
used to mask the individual modulated, normalized GM 
images and extract the average volume within each ROI 
using the REX toolbox (http://web.mit.edu/swg/software. 
htm). A linear regression analysis was performed at 
every ROI to remove the effects of age and total brain 
volume. The residuals of this regression were then sub- 
stituted for the raw ROI volume values [22,27,36]. 

Construction of structural correlation network 

The extracted residual volumes of all 90 anatomical 
ROIs were used for construction of structural correlation 
networks. For each group, a 90 x 90 association matrix R 
was generated with each entry r t j defined as the Pearson 
correlation coefficient between residual volumes of 
regions i and across subjects. From each association 
matrix, a binary adjacency matrix A was derived where 
a t j was considered 1 if r t j was greater than a specific 
threshold and zero otherwise. The diagonal elements of 
the constructed association matrix were also set to zero. 
The resultant adjacency matrix represented a binary un- 
directed graph G in which regions i and ; were con- 
nected if gij was unity. Therefore, a graph was 
constructed with N=90 nodes (anatomical ROIs), with a 
network degree of E equal to number of edges (links), 
and a network density (cost) of D = E/[(N x (N-l))/2] 
representing the fraction of present connections to all 
possible connections. Since thresholding the association 
matrices of different groups at an absolute threshold 
results in networks with a different number of nodes 
(and degrees) that might influence the network measures 
and reduce interpretation of between group results [37], 
two approaches were implemented for thresholding the 
constructed association matrices based on previous stud- 
ies [26-28]: (1) Thresholding the constructed association 
matrices at a minimum network density in which all 
nodes become fully connected in the brain networks of 
both groups; (2) Thresholding the constructed associ- 
ation matrices at a range of network densities and com- 
paring the network topologies across that range. 

Network analysis 

Small-worldness The small-worldness of a complex 
network, as described above, has two key metrics: the 
clustering coefficient C and the characteristic path 
length L of the network. The clustering coefficient of a 
node is a measure of the number of edges that exist be- 
tween its nearest neighbors. The clustering coefficient of 
a network is thus the average of clustering coefficients 
across nodes and is a measure of network segregation. 
The characteristic path length of a network is the aver- 
age shortest path length between all pairs of nodes in 
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the network and is the most commonly used measure of 
network integration [38]. To evaluate the topology of 
the brain network, these parameters must be compared 
to the corresponding mean values of a random graph 
with the same number of nodes, total edges, and degree 
distribution as the network of interest [39,40]. Thus, we 
obtained the small-worldness index of a network as 
[C/C rand ]/[L/L rand ] where C rand and L rand are the mean 
clustering coefficient and the characteristic path length 
of the random network [20]. In a small -wo rid network, 
the clustering coefficient is significantly higher than that 
of random networks {C/C rand ratio greater than 1) while 
the characteristic path length is comparable to random 
networks (L/L rand ratio close to 1). 

Regional network measures We also investigated the 
nodal characteristics of the constructed structural net- 
works to identify differences in regional network mea- 
sures between groups. Nodal betweenness centrality and 
nodal degree were calculated for each of the anatomical 
ROIs for the networks thresholded at minimum density 
with full-connectivity. Nodal betweenness centrality is 
defined as the fraction of all shortest paths in the net- 
work that pass through a given node and is used to de- 
tect important anatomical or functional connections. 
Nodes that bridge disparate parts of the network have a 
high betweenness centrality [38]. On the other hand, 
nodal degree is defined as the number of connections 
that a node has with the rest of the network and is con- 
sidered a measure of interaction of a node, structurally 
or functionally, with the network. 

Network hubs Hubs are crucial components for effi- 
cient communication in a network. Hubs are not only 
considered as important regulators of information flow 
but also play a key role in network resilience to insult 
[38]. We considered a node as a hub if its betweenness 
centrality was at least 2SD higher than mean network 
centrality [26]. 

Comparing correlation strengths between groups To 

test the significance of differences in overall inter- 
regional correlation of gray matter volume between net- 
works, the correlation coefficients were converted to z 
values using Fishers r-to-z transform. This transform- 
ation resulted in values that were approximately normally 
distributed [28]. Then, a two sample £-test was used to 
test the significance of the difference in mean overall cor- 
relation between groups. 

Comparing network measures between groups In 

order to test the statistical significance of the between- 
group differences in network topology and regional net- 
work measures, a non-parametric permutation test with 



1000 repetitions was used [26,28]. In each repetition, 
the calculated residual volumes of each participant were 
randomly reassigned to one of the two groups so that 
each randomized group had the same number of sub- 
jects as the original groups. Then, an association matrix 
was obtained for each randomized group. The binary 
adjacency matrices were then estimated by thresholding 
the association matrices at a range of network densities. 
The network measures were then calculated for all 
the networks at each density. The differences in net- 
work measures between randomized groups were 
then calculated resulting in a permutation distribu- 
tion of difference under the null hypothesis. The ac- 
tual between-group difference in network measures 
was then placed in the corresponding permutation 
distribution and a two-tailed p-value was calculated 
based on its percentile position [27]. 

We used the Brain Connectivity Toolbox [38] for 
quantification of network measures in addition to in- 
house software, Graph Analysis Toolbox (http://nnl. 
stanford.edu/Tools.html), for comparing the structural 
networks. Brain Net Viewer (http://www.nitrc.org/ 
projects/bnv/) was used for visualization of the 
graphs. 

It should be noted that we could not statistically ex- 
plore the relationship between cognitive-behavioral mea- 
sures and the extracted network properties. It is because 
the extracted network measures are group-specific and 
can not be quantified for individuals. Therefore, we 
could not perform statistical analysis (e.g. correlation 
analysis) to relate network properties and cognitive- 
behavioral outcomes. 

Results 

Comparing networks at minimum density 
with full-connectivity 

The minimum network density in which all nodes be- 
came fully connected in the structural networks of both 
groups was 0.184. The association matrices for each 
group and corresponding binary adjacency matrices 
thresholded at the minimum density are shown in 
Figure 1. The correlation network of the BC group 
showed lower overall correlation strength than the CON 
group (t = 5.9, P <0.001). Topological measures for the 
networks of both groups thresholded at the density of 
0.184 are given in Table 2. In both groups, the normal- 
ized characteristic path length of the network was close 
to 1 and the normalized clustering coefficient was 
greater than 1 which results in a small-worldness par- 
ameter of greater than 1. Thus, consistent with previous 
studies, the gray-matter correlation network of both 
groups followed a small-world property. While no sig- 
nificant difference in characteristic path length was 
observed between networks (P = 0.22), the normalized 
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Figure 1 Association and binary adjacency matrices; association matrices for (A) BC and (B) CON groups, binary adjacency matrices for 

(C) BC and (D) CON groups (connected regions are shown in red). The association matrices show the highest connectivity between regions 
of interest as well as for inter-hemispheric regions. 



clustering coefficient was significantly higher in CON 
network (P<0.05). Consequently, the small-worldness 
parameter was found to be greater in the CON network 
than in the BC network and the observed difference was 
marginally significant (P = 0.08). 



Table 2 Network measures at minimum density of 0.184 
with full connectivity 





BC 


CON 


p-value 


Density 


0.184 


0.184 




Mean clustering coefficient 


0.485 


0.517 


0.16 


Characteristic path length 


2.11 


2.2 


0.23 


Normalized clustering 


2.65 


3.06 


0.03 


Normalized path length 


1.15 


1.20 


0.22 


Small-worldness 


2.30 


2.54 


0.08 



Comparing networks across a range of density 

In order to investigate changes in the network topology 
as a function of network density (cost), we thresholded 
the constructed association matrices at a range of net- 
work densities. Changes in network parameters are 
depicted in Figure 2. The figures show that both net- 
works follow a small-world organization across the 
range of densities. The results of comparison of the 
two networks as a function of network density are 
shown in Figure 3. While there was no significant dif- 
ference in network characteristic path lengths between 
groups (same level of integration), the clustering coef- 
ficient was significantly higher (more segregated net- 
work) in the CON network in various densities 
resulting in higher small-worldness compared to the 
BC group across a range of densities. 
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Figure 2 Global network measures; A) clustering coefficient, B) 
characteristic path length and C) small-worldness of BC and 
CON networks. The figures show that both networks follow a 
small-world organization across the range of densities; the characteristic 
path length is close to 1 (for the fully connected network, i.e. 
density > 0.19) and the clustering is greater than 1 in different densities. 



Below the minimum density of full connectivity 
(0.184), the BC network fragments sooner that the CON 
network and thus the normalized path length of the BC 
network increases relatively faster compared with that of 
CON network. This results in a relatively large drop in 
the difference in normalized path length between groups 
in densities below 0.184. 

Regional network measures 

As shown in Figure 4, the CON group demonstrated 
higher nodal betweeness/degree in several frontal and 
temporal regions. The BC group demonstrated higher 
nodal betweeness/degree primarily in parieto-occipital 
regions but also in the anterior cingulate and left 
thalamus. 

Network hubs 

In the BC group, network hubs were identified at left an- 
terior cingulate, right inferior parietal lobule, and right 
supramarginal gyrus. In CON group, hubs were found in 
right middle frontal gyrus, left superior frontal gyrus, 
right insula, right precuneus, and left middle temporal 
gyrus (Figure 5). 

Discussion 

We applied graph theoretical analyses to compare gray 
matter volume correlation networks of female primary 
breast cancer (BC) survivors treated with chemotherapy 
and female healthy controls (CON). We found an alter- 
ation in small-world characteristics of the brain struc- 
tural network in the BC group; an observation that 
confirms our hypothesis suggesting changes in large- 
scale brain network properties in BC survivors treated 
with chemotherapy. To our knowledge, this is the first 
report of altered large-scale network properties asso- 
ciated with BC and chemotherapy. 

Between-group differences in network topology 

While the structural correlation networks of both groups 
followed a small-world organization, the results revealed 
a decrease in the small-world index of the BC group 
compared to controls. The characteristic path lengths of 
the BC and CON group networks were not different 
whereas the clustering coefficient was significantly smal- 
ler in the BC network than in the CON network. 
Decreased clustering coefficient and small-world index 
in the BC group was observed not only at the minimum 
network density but also over a range of density 
thresholds. 

These findings suggest that the structural correlation 
networks of BC patients tend to have a more rando- 
mized configuration compared to CON group associated 
with weaker regional connectivity and disrupted global 
organization [38]. Given that BC tends to be a disease 
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associated with aging, altered network topology could 
make women with BC more vulnerable to age-related 
cognitive decline. Recent studies suggest that older BC 
patients show increased chemotherapy-related cognitive 
deficits [4,5]. Additionally, this same profile of network 



Figure 3 Between-group differences in global network 
measures; between-group differences in A) clustering 
coefficient B) characteristic path length and C) small-worldness 
parameter. The vertical arrows represent the densities where the 
difference is statistically significant at P<0.05. It shows that while 
there is no significant difference in network characteristic path 
lengths between groups (same level of integration), the clustering 
coefficient is significantly higher (more segregated network) in the 
CON network in various densities resulting in higher small-worldness 
compared to the BC group across a range of densities. Squares 
represent control minus BC group measure, dashed lines show the 
95% confidence interval and the dotted line shows the mean for the 
random graph distribution. 



topology (i.e. decreased clustering and small-world 
index) has been associated with Alzheimer's dementia 
[41], The interface between aging, breast cancer and 
cognitive decline will be a critical area of continued 
investigation. 

The observed lack of densely interconnected neighbors 
in the network of BC group is also supported by the 
lower overall correlation strength in the association 
matrix of BC group compared to CON group. Together, 
the results of network-level analysis corroborates previ- 
ous structural neuroimaging findings that have demon- 
strated a diffuse pattern of atrophy in gray matter 
volumes in BC patients [15,16] and further suggests an 
alteration in the coordinated patterns of brain morph- 
ology in the structural network of these patients. 

One potential mechanism underlying this network- 
level alteration is white matter damage. There is a large 
body of evidence suggesting that global gray matter atrophy 
is associated with focal and global white matter damage 
due in part to the transection of axons and subsequent 
retrograde neuronal loss [42,43]. Animal studies show that 
chemotherapy suppresses neural progenitor cell prolifera- 
tion responsible for white-matter tract integrity and 
cortico-cortical connections [44-48]. Thus, the observed 
network-level alteration in structural correlation network 
of BC patients might arise from neurotoxic effects of 
chemotherapy on cortico-cortical connections. This idea 
is also supported by a recent diffusion weighted imaging 
study involving chemotherapy-treated BC survivors that 
reported a diffuse pattern of microstructural white matter 
damage [14], 



Between-group differences in regional network measures 

High degree/betweenness nodes were observed mainly 
in frontotemporal regions in the CON network and in 
parietal, cingulate, and occipital regions in the BC net- 
work. Nodes with high degree/centrality in the structural 
network identify regions that are highly interactive and 
have the potential to participate in a large number of 
functional interactions [49]. The results indicate that the 
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Figure 4 Difference in regional network characteristics relative to random networks; A) Controls showed greater network betweenness 
centrality in left SFG, left PrCG, right PCUN, right REC and right ITG while the BC group showed greater betweenness in right ANG, 
right CALC, left ACQ right IPL, and left THL. B) Controls demonstrated higher network degree in left MFG, left MFOr, bilateral SFG, left 
PrCG, right ITG and right MTG. The BC group showed greater degree in bilateral ANG, left IOG, right IPL and left SMG. Squares represent control 
minus BC group measure, dashed lines show the 95% confidence interval and the dotted line shows the mean for the random graph distribution. 
Regions that showed significantly higher/lower degree in CON relative to BC are shown in pink/cian color on the ICBM152 brain template. 
Abbreviations are used as follow: L: left hemisphere; AMYG: amygdala; ANG: angular gyrus; CALC: calcarine fissure; CN: caudate nucleus; ACC: 
anterior cingulate; MCC: mid-cingulate; PCC: posterior cingulate; CUN: cuneus; IFOp: inferior frontal gyrus, opercular part; IFOr: inferior frontal gyrus, 
orbital part; IFTr: inferior frontal gyrus, triangular part; MedFOr: medial frontal gyrus, orbital part; MFG: middle frontal gyrus; MFOr: middle frontal 
gyrus, orbital part; SFG: superior frontal gyrus; MedSF: superior frontal gyrus, medial part; SFOr: superior frontal gyrus, orbital part; FG: fusiform 
gyrus; HSHL: heschl gyrus; HIPP: hippocampus; INS: insula; LNG: lingual gyrus; IOG: inferior occipital gyrus; MOG: middle occipital gyrus; SOG: 
superior occipital gyrus; OFB: olfactory cortex; PLD: lenticular nucleus, pallidum; PCL: paracentral lobule; PHIP: parahippocampal gyrus; IPL: inferior 
parietal lobule; SPL: superior parietal lobule; PoCG: postcentral gyrus; PrCG: precentral gyrus; PCUN: precuneus; PUT: putamen; REC: gyrus rectus; 
RLN: rolandic operculum; SMA: supplementary motor area; SMG: supramarginal gyrus; ITG: inferior temporal gyrus; MTG: middle temporal gyrus; 
MTP: middle temporal pole; STP: superior temporal pole; STG: superior temporal gyrus; THL: thalamus. 



network of BC group has less highly-interactive nodes in 
the frontotemporal regions. 

Reduction in gray and white matter volumes in fronto- 
temporal regions of patients with BC before and after 
treatment has been shown in a recent study [15]. A 
number of functional neuroimaging studies have also 



reported increased activation in the frontotemporal net- 
work in these patients suggesting necessity of recruiting 
a compensatory neural mechanism to complete certain 
tasks, or decreased efficiency of these systems [10,50,51]. 
The observed changes in the centrality of the frontotem- 
poral nodes in the network of BC patients is consistent 
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with these previous results and may explain the impair- 
ment of these patients in memory and executive 
functioning. 

There are several mechanisms that could explain the 
observed change in the distribution of highly influential 
nodes in BC group with less highly- interactive nodes in 
the prefrontal and temporal cortices. Certain genetic and 
immunologic factors have been proposed as important 
contributors to cognitive decline in patients with BC [3]. 
For example, apolipoprotein E (APOE) is a complex pro- 
tein that has an important role in neuronal repair and 
plasticity after injury [52], The relationship between this 
genotype and adverse cognitive performance and reduce 
mesial temporal (hippocampal) volume has been shown 
in long-term survivors of BC [53]. Also, cytokines that 
become dysregulated during cancer and/or chemother- 
apy treatment have been shown to affect cognitive func- 
tion in healthy adults [54]. Cytokines have important 
roles in central nervous system functioning including 
modulation of dopamine and other neurotransmitters 
[55] that affect the frontotemporal network [56]. The 
effects of genotypic variations and cytokine regulation 
on large-scale networks have not been established and 
require further investigation. 

The BC group had significantly more post-menopausal 
women than the healthy comparison group. This is 



expected given that chemotherapy frequently causes 
early menopause [57] and tamoxifen can increase meno- 
pausal symptoms [19], both by blocking or reducing es- 
trogen. Estrogen has a marked effect on brain function, 
particularly in frontotemporal regions and their con- 
comitant executive and memory functions [58]. Previous 
studies suggest that tamoxifen, in addition to chemo- 
therapy, can have detrimental effects on cognitive and 
neurobiologic outcomes in BC survivors [59,60]. The 
present findings may indicate potential secondary neu- 
roendocrine effects on large-scale brain networks. 

Network hub analysis 

There was a difference in the number and distribution 
of hubs between the BC and CON groups. We found six 
network hubs in the CON group mainly across frontal 
and temporal regions compared to three hubs in the BC 
group primarily in parietal and cingulate areas. The 
identified hubs in the CON group are consistent with 
the results of previous graph-theoretical analysis involv- 
ing healthy adults [26,28]. The present findings indicate 
that the network of BC patients does not show expected 
prefrontal and temporal hubs suggesting network altera- 
tions involving these regions, again which are critical for 
the executive function and memory skills commonly dis- 
rupted in BC. 



Hosseini et al. BMC Neurology 2012, 12:28 
http://www.biomedcentral.com/1471-2377/12/28 



Page 10 of 11 



Additionally, the number of identified hubs in the BC 
network was half the number of hubs in the CON net- 
work. Since hubs are crucial components for efficient 
communication across the network, this observation 
suggests a less efficient information transfer in the net- 
work of BC patients. These results are consistent with 
our previous findings showing that BC survivors can 
often perform similarly to controls on cognitive tasks 
but require significantly more neural resources to do so 
[50]. Standardized neurocognitive tests may lack sensi- 
tivity to these subtle changes [2] as most were designed 
to assess specific cognitive skills. Based on the present 
findings, measures of reaction time or parallel processing 
should be included in evaluations of BC-related cogni- 
tive outcome. Quantitative process measures might also 
be more effective as these assessments examine the effi- 
ciency of the patient s performance strategies in addition 
to performance accuracy [61]. 

This study has several limitations. First, this was a 
cross-sectional study and therefore, we could not dir- 
ectly test the dissociated effects of cancer and chemo- 
therapy on network measures. Future studies can 
address this issue through longitudinal evaluations of 
network measures in patients with BC. In addition, fu- 
ture studies might investigate whether the alteration in 
gray matter network is due to any cancer (and not just 
breast cancer). It is also important to dissociate between 
cancer-related physiological changes and cancer effects 
on mental states. Second, the network measures were 
identified by calculating the correlations of gray matter 
volume across subjects. While this methodology is com- 
mon for investigating between-group differences in brain 
structural network, it lacks the ability to explore the in- 
dividual differences in network parameters. Specifically, 
we could not explore the effect of treatment period, regi- 
men, age or other individual differences on network 
parameters or investigate the effects of altered network 
parameters on cognitive-behavioral outcome. 

Conclusions 

Despite these limitations, these findings increase our 
understanding of BC and chemotherapy-related cogni- 
tive impairment by demonstrating alterations in specific 
network properties. Survivors of BC treated with chemo- 
therapy showed decreased regional connectivity (cluster- 
ing coefficient) and global network organization (small- 
worldness) and integration (hubs) suggesting reduced 
robustness and efficiency of the network. A pattern of 
frontotemporal abnormality was also noted. These 
results contribute novel insights regarding the neurobio- 
logy mechanisms underlying cognitive deficits in these 
patients and highlight critical areas for future research 
related to the potential vulnerability of BC patients to 



the neurologic effects of aging and other forms of struc- 
tural damage. 
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